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In this communication, we report results of three-dimensional hydrodynamic compu- 
q ■ tations, by using equations of state with a critical end point as suggested by the lattice 
QCD. Some of the results are an increase of the multiplicity in the mid-rapidity region 
P-i! and a larger elliptic-flow parameter w 2 . We discuss also the effcts of the initial-condition 
Q_i' fluctuations and the continuous emission. 

1. INTRODUCTION 

- T— I | 

^ ■ Nowadays, it is widely accepted that hydrodynamics is a successful approach for de- 
scribing the bulk of the collective flow in high-energy nuclear collisions [ I]. The basic 
assumption in hydrodynamical models is the local thermal equilibrium. Once this con- 
dition is satisfied, all the thermodynamical relations should be valid in each space-time 
point. The properties of the matter formed in high-energy collisions are then specified by 
some equations of state (EoS). Thus, one of the main objects of hydrodynamical approach 
is to determine which are the EoS that consistently reproduce the observed quantities. 

In high-energy nucleus-nucleus collisions, one often uses EoS with a first-order phase 
transition, connecting a high-temperature QGP phase of the system with a low-temperature 
hadronic phase. However, lattice QCD studies showed that the transition line between 
QGP and hadron phase has a critical end point and for small net baryon surplus the 
transition is of crossover type [E]. So, we would like to learn what are the consequences 
of these results on the hydrodynamics and on the observable quantities. 

We shall begin showing, in the next Section, how the critical end point could be imple- 
mented for the sake of phenomenological computations. Besides the EoS, the ingredients 
of any hydrodynamic approach are the equations of motion, the initial conditions and 
some decoupling prescription. In Sec. EH we shall discuss how the initial conditions are 
chosen in our studies; then, how we solve the hydrodynamic equations; and finally which 
is the decoupling prescription. In Sec. we shall show some of the results of our studies. 
Finally, summaries of conclusions and outlook are given in Sec. |S] 
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2. PARAMETRIZATION OF EQUATIONS OF STATE 



As mentioned above, one often introduces EoS with a first-order phase transition, con- 
necting a QGP phase of the system, usually described by the MIT bag model, with a 
hadronic phase, depicted as a resonance gas. A detailed account of such EoS may be 
found, for instance, in [|2|. We start from these EoS, in order to get a phenomenological 
parametrization of those suggested by lattice QCD. 

Let us denote by Pq the pressure given by the MIT bag model and Ph the one cor- 
responding to the hadronic resonance gas. Given a value //& of the baryonic chemical 
potential, we write for the pressure P the equation (see Figure HJ) 



(P - P Q ){P - P H ) = 6(ji b ) 
where 

S(fJ>b) = <5 exp -(Hb/t*c) 



with fi c = critical chemical potential. 



By solving the equation above and using thermodynamical relations, we obtain 

9 c 

P = \P H + (1-\)P Q + 

s = Xs H + (1 - A)s Q , 
n b = Xn H + (1 - \)n Q - 

e = Xe H + (l- X)€q 



(P Q -P H ) 2 + 45 
2 (fi/fil) 5 



2 [1 + ( M / Mc 



2 + 45 
) 2 1 6 



(P Q - P H f + 45 



where A = - 



1 - (P Q - P H )/y/(PQ-P H ) 2 + 4S 



(1) 
(2) 

(3) 

(4) 
(5) 

(6) 
(7) 



Observe that if 5q = 0, we recover the EoS with the first-order phase transition described 
above. In the discussion below, we call them 10PT EoS. As seen in Figure^ when <5(/if>) 7^ 
0, the transition from hadron phase to QGP is smooth. In the right-hand side of Eq. ([T|) . we 

could choose some function which be- 
comes exactly for fi b > fi c to guarantee 
the first-order phase transition there, but 
for practical purpose this is not necessary. 
As will be shown below, our choice repre- 
sented by Eq.(J2J) is enough. We shall des- 
ignate the equations of state given above, 
with 5 ^ 0, CP EoS. 

Let us compare in Figure 121 below, the 
temperature dependences of the energy 
and entropy densities, e and s, and the 
pressure P, given by the two sets of EoS 
defined above. One can see that the cross- 




Figure 1. Illustration of Eqs.((TJ) and (0). 
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Figure 2. A comparison of e(T), s(T) and P(T) as given by our parametrization with a 
critical point (solid lines) and those with a first-order phase transition (dashed lines). 



M=0 GeV 



u=0.2 GeV 



u=0.4 GeV 




Figure 3. Plots of s/e and P as function of e for the two EoS shown in Figure El 
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over behavior is correctly reproduced by our parametrization for CP EoS, while finite 
jumps in s(T) and s(T) are exhibited by 10PT EoS, when T crosses the transiton tem- 
perature. It is also seen, as mentioned above, that at ~ 0.4 GeV the two EoS are 
indistinguishable. 

Now, since in a real collision what is directly given is the energy distribution at a 
certain initial time (besides baryon number distribution, charge distribution, strangeness 
distribution, etc.), whereas the temperature is defined with the use of the former, it would 
be nice to compare the two sets of EoS, by plotting several quantities as function of e. 
We do this in Figure H3 One immediately sees there some remarkable differences between 
the two sets of EoS: i ) naturally the pressure is not constant for CP EoS in the crossover 
region; ii ) moreover, the entropy is larger. We will see in Sec. 4 that these characteritics 
affect the observed quantities in non-negligible way. 



3. OTHER INGREDIENTS OF OUR HYDRODYNAMIC MODEL 

Besides the equations of state, the other ingredients of a hydrodynamic model are the 
initial conditions, the equations of motion and some decoupling prescription. Here we 
shall discuss how these elements are chosen in our studies. 

3.1. Initial Conditions 

In usual hydrodynamic approach, one assumes some highly symmetric and smooth 
initial conditions (IC). However, since our systems are small, large event-by-event fluctu- 
ations are expected in real collisions, so this effect should be taken into account. Remark 
that this might happen even if the impact parameter could be maintained fixed. 

Many simulators, based on microscopic models, e.g. HIJING [H], VNI [E|, URASiMA 
[El, NeXuS [|7j, • • •, show such event-by-event fluctuations. As an example we show in 
Figure E] the energy density for central Au+Au collisions at 130A GeV, given by NeXuS 
simulator [|7j, at mid-rapidity. In our approach, we use both fluctuating and averaged IC. 
Some consequences of such fluctuations have been discussed elsewhere. We shall discuss 
some others in Sec HI 




Figure 4. The energy density is plotted in units of GeV/fm 3 at the initial instant. Left: 
one random event. Right: average over 30 random events (corresponding to the smooth 
initial conditions in the usual hydro approach). 
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3.2. Equations of Motion 

The equations of motion of hydrodynamics are the continuity equations expressing the 
energy-momentum conservation, the baryon-number conservation, and other conservation 
laws, corresponding to several charges. Here, for the sake of simplicity, we shall consider 
only the energy-momentum and the baryon number. Since our IC are entirely arbitrary, 
without any symmetry, as discussed above, the only way to solve the equations is through 
numerical computations. We have developed a numerical code called SPheRIO (Smoothed 
Particle hydrodynamic evolution of Relativistic heavy IOn collisions) [|Hj, based on the 
so called Smoothed- Paricle Hydrodynamics (SPH) algorithm [|5J. 

3.3. Decoupling Prescription 

In hydrodynamic treatment of high-energy nuclear collisions, one often assumes decou- 
pling on a sharply defined hypersurface, usually characterized by a constant temperature 
Tf Q . We call this Sudden Freeze Out (FO). However, our systems are small, so particles 
may escape from a layer with thickness comparable with the systems sizes. We have 
proposed an alternative decoupling prescription that we call Continuous Emission (CE) [ 
ITU] which, as compared to the usual sudden freeze out, we believe closer to what happens 
in the actual collisions. We introduce, at each space-time point x^, a certain momentum- 
dependent escaping probability 



V(x, k) = exp 



/oo 
p(x') ov dr' 



(8) 



To implement this prescription in our SPheRIO code, we had to introduce some ap- 
proximation to make the computation practicable. First, we take V on the average, i.e., 

V(x,k) -> (V(x,k)) = V(x) (9) 

and then we approximate linearly the density p(x') = a s(x') in Eq.fjHJ). Thus, 



V(x, k) - V{x) = exp (-k J^j^J > ( 10 ) 

where k = 0.5 a (av) is estimated to be 0.3 , corresponding to (av) ~ 2 fm 2 . 

We show, in Figure the time (r) evolution of the probability V, estimated by this 
expression, in the mid-rapidity plane for the most central Au+Au collisions at 130 A GeV. 
For comparison, we show, in Figure® the corresponding time evolution of the temperature 
T. One sees that, whereas in Sudden Freeze Out, particles are emitted from a constant T 
line of Figure |H1 in Continuous Emission, they are emitted according to the probability 
V, so from a difuse space region and during a larger time interval. It will be shown in 
Sec. HJthat this difference gives important changes in some observables. 
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x(fm) 

Figure 5. Time evolution of the probability 
as given by Eq.()10|) for the most central 
Au+Au collisions at 130 A GeV, in the mid- 
rapidity plane for averaged IC. 




x(fm) 

Figure 6. Corresponding time evolution of 
the temperature. 



4. RESULTS 



Let us now show results of computation of some observables, as described above, 
for Au+Au collisions at 200A GeV. We start computing the pseudo-rapidity and the 
transverse-momentum distributions for charged particles, to fix the parameters. Then, 
the elliptic-flow parameter t> 2 and HBT radii are computed in fit-parameter free way. 

4.1. Pseudo-rapidity distribution 

In Sec. 121 we showed that the inclusion of a critical end poit in the first-order phase- 
transition line increases the entropy per energy, as compared with 10PT EoS. This means 
that, given the same total energy, the multiplicity is larger for CP EoS case than for 10PT 
EoS one. Figure [7| shows clearly that this happens, especially in the mid-rapidity region. 

Now, once the equations of state are chosen, what is the effect of the fluctuating initial 
conditions for the same decoupling prescription? This effect has already been discussed 
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Figure 7. rj distributions for the most central 
Au+Au at 200A GeV. A comparison of CP 
EoS (solid line) vs. 10PT EoS (dashed line). 
The data are from PHOBOS Collab.[Q2]. 



Figure 8. rj distributions for the most cen- 
tral Au+Au at 200A GeV, computed with 
CP EoS, with three different combinations 
of IC and decoupling prescriptions. 
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earlier and shown in Figure |H1 namely, the multiplicity is smaller if the IC fluctuation 
is taken into account and average computed after the decoupling, as compared with the 
result obtained with smooth averaged IC. The curve with continuous emission is also 
shown there, reproducing equally well the data. Here, the parameter k has been fixed as 
explained in Subsec. 



4.2. Transverse-Momentum Distribution 

As discussed in Sec. 121 since the pressure does not remain constant in the crossover 
region, we expect that the transverse acceleration is larger for CP EoS, as compared with 
10PT EoS case. In effect, Figure |H1 does show that px distribution is flatter for CP EoS, 
but the difference is small. 

We show, in Figure ITUl three different combinations of IC and decoupling prescriptions, 
corresponding to CP EoS. The curve with the event-by-event fluctuating IC is flatter than 
the one corresponding to the averaged IC, both with sudden freezeout, probably because 
the initial expansion in the former is more violent, due to the bumpy structure with high- 
density blobs, as seen in Figure |U in this case. The freezeout temperature suggested by 77 
and pt distributions turned out to be Tf ~ 135 — 140 MeV. 
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Figure 9. px distributions for the most cen- 
tral Au+Au at 200A GeV. A comparison of 
CP EoS (solid line) vs. 10PT EoS (dashed 
line). The data are from PHOBOS Collab.[ 



Figure 10. pt distributions for the most 
central Au+Au at 200A GeV, computed 
with CP EoS, with three different com- 
binations of IC and decoupling prescrip- 
tions. 



4.3. Elliptic-Flow Parameter v 2 

Here, we show our results for the pseudo-rapidity distribution of the elliptic-flow pa- 
rameter V2 for Au+Au collisions at 200A GeV. As seen in Figure CP EoS gives larger 
vi , as a consequence of larger acceleration in this case as discussed in Sec|21 Notice that 
the continuous emission makes the curves narrower, as a consequence of earlier emission 
of particles, so smaller acceleration, at large- 1 77 1 regions. See more detils in[E]- When the 
IC fluctuations are taken into account, the resulting fluctuations of t>2 become large, as 
seen in Figures and It would be nice to measure such a v 2 distribution, which would 
discriminate among several microscopic models for the initial stage of nuclear collisions. 
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Figure 11. Pseudo-rapidity distribution of 
t> 2 for charged particles in the centrality 
(15 - 25)% Au+Au at 200A GeV, computed 
with event-by-event fluctuating IC. The ver- 
tical bars indicate dispersions. The data are 
from PHOBOS Collab.[E|. 



Figure 12. Event-by-event v% distribution 
in the interval 0.48 < i] < 0.95 and corre- 
sponding to the same events as in Figure 

HU 



4.4. HBT Radii 

Here, we show our results for the HBT radii, in Gaussian approximation as used in 
experimental data analyses, for the most central Au+Au collisions at 200A GeV. As seen 
in Figures Unillll and [El the differences between CP EoS results and those for 10PT EoS 




Figure 13. h? depenence of HBT radius Rl 
for pions in the most central Au+Au at 200A 
GeV, computed with event-by-event fluctu- 
ating IC. The data are from PHENIX Col- 
lab.[H5|. 



Figure 14. kx depenence of HBT radius 
R s for pions in the most central Au+Au 
at 200A GeV, computed with event-by- 
event fluctuating IC. The data are from 
PHENIX CoUab.[Hnj. 
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Figure 15. kx depenence of HBT radius R Q 
for pions in the most central Au+Au at 200A 
GeV, computed with event-by-event fluctu- 
ating IC. The data are from PHENIX Col- 
lab, [ng. 



Figure 16. depenence of the ra- 

tio R Q /R S for pions in the most central 
Au+Au at 200A GeV, computed with 
event-by-event fluctuating IC. The data 
are from PHENIX Collab.[E]. 



are small. This was expected for the case of Rl ■ For R s , and especially for R , one sees 
that CP EoS combined with continuous emission gives steeper kx dependence, closer to 
the data. However, there is still numerical discrepancy in this case. 

5. CONCLUSIONS AND OUTLOOKS 

In this work, we have introduced a phenomenological parametrization of lattice-QCD- 
inspired equations of state, which presents a first-order phase transition at large baryonic 
chemical potential and a crossover behavior at smaller chemical potential. By using the 
initial conditions generated by NeXuS event simulator and SPheRIO code for solving the 
hydrodynamic equations, some observables were computed and studied the effects of such 
EoS. Some of the conclusions are: 

1. The multiplicity becomes larger for these equations of state in the mid-rapidity; 

2. The pt distribution becomes flatter. However, the differece is small; 

3. Larger t> 2 • Continuous Emission makes the rj distribution narrower; 

4. HBT radii slightly closer to data. 

In our calculations, the effect of the continuous emission on the interacting component 
has not been taken into account. A more realistic treatment of this effect probably 
makes R smaller, because the duration for particle emission becomes smaller in this case. 
Another improvement we should make is the approximations we used for Eq.fjHJ). Cascade 
treatment of this part is probably a better alternative. 
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